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Abstract. We report the measurement of the fluctuations of a transmon qubit 
through the noise spectrum of the microwave signal that measures its state. The 
amplitude of the Lorentzian noise power spectrum allows to determine the average 
qubit excitation, in agreement with the estimated thermal radiation reaching the 
sample. Its width yields the qubit energy relaxation rate which decreases with 
temperature, contrary to the predictions for a two-level system solely coupled to 
thermal radiation. This indicates the existence of another non-radiative energy 
relaxation channel for the qubit. 
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Superconducting qubits pj are promising candidates for implementing a solid- 
state quantum processor. Over the last years, substantial improvements have been 
made in the coherence times [2j [3], fidelity of single-qubit gates [4j, readout procedures 
[SI [6] and entanglement of several qubits [7]. In a recent experiment a simple 
quantum algorithm was operated on a two-qubit elementary processor [8]. One of 
the requirements for the implementation of larger scale quantum algorithms [9j is that 
the qubit registers should be properly initialized at the beginning of each computation, 
with all qubits lying in their ground state. In most superconducting qubit experiments, 
the initialization is simply realized by waiting long enough before each experimental 
sequence for the system to reach thermal equilibrium. At cryogenic temperatures in 
the 10 — 30 mK range and for typical qubit resonance frequencies of a few GHz, there 
is indeed at thermal equilibrium a small (typically less than 1%) probability of finding 
the qubits in the excited state, which is usually considered negligible. However, given 
the recent improvement of the overall fidelity of single- and two-qubit gates, the effect 
of even small thermal fluctuations will require to be considered more quantitatively 
in the near future. Moreover it is well known in mesoscopic physics that the effective 
temperature of an electrical degree of freedom such as a superconducting qubit can 
be in some cases much larger than the temperature of the cryostat, because it can 
be strongly coupled to out-of-equilibrium electromagnetic radiation coming from the 
measuring leads, while only weakly to the phonon bath. It is thus important to be 
able to measure precisely the average excited state population of a single qubit. 

Here we propose and demonstrate a method to determine this thermal excited 
state population in a circuit Quantum Electrodynamics (cQED) setup [10], where a 
Cooper-pair box qubit of the transmon type [HI [3] is coupled to a coplanar waveguide 
resonator (CPW). The two qubit states shift differently the resonator frequency, so 
that the phase of a microwave signal reflected by the resonator allows a non-destructive 
readout as demonstrated in numerous experiments [12] • In the present setup as in 
most cQED experiments, it is not possible to readout the qubit state in one single 
experimental sequence due to insufficient signal-to-noise ratio. Note however that 
such a single-shot readout was recently obtained in cQED by using a non-linear CPW 
resonator [5]. The usual method for reading out the qubit state, using ensemble- 
averaged measurements of the microwave signal, does not directly provide an absolute 
measurement of the qubit excitation. However, thermal fluctuations of the qubit state 
are responsible for a measurable phase noise in the microwave signal reflected by the 
resonator, with a characteristic Lorentzian power spectrum. In this work we report 
the observation of this thermal noise and we use it to determine the effective qubit 
temperature. We note that a related measurement was performed on an ensemble of 
nuclear spins measured by a SQUID amplifier [13]. 

The complete experimental setup is shown in Fig [TJ The transmon has its 
two lowest energy eigenstates \g) and |e) separated by u} ge /2ir =5.304 GHz. It is 
capacitively coupled with strength g/2w = 45±2 MHz to a superconducting resonator 
with resonance frequency uj c /2ir =5.796 GHz and bandwidth BW =30.3 MHz, which 
serves as qubit readout. With these parameters, the qubit is sufficiently detuned 
from the resonator for their interaction to be well described by the dispersive 
Hamiltonian H — hxna z , where x ls the dispersive coupling constant and n is the 
intra-resonator photon number operator. The resonator frequency is thus shifted by 
±x/27r = 1.75 MHz when the qubit is in \g) or \e) respectively. A continuous microwave 
tone of frequency u c /2ir sent to the resonator input from source V m acquires a qubit 
state-dependent phase shift which allows a continuous and non-destructive monitoring 
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Figure 1. Microwave setup used for measurements. A microwave signal V m 
is sent through the input line containing several attenuators and filters at each 
temperature stage to the input port of the resonator (shown in green). The 
reflected signal is separated from the input one by a circulator and goes through 
a filter and two isolators before reaching the cryogenic amplifier (gain 37 dB and 
noise temperature 2.6 K). The signal is then demodulated at room temperature 
with a homodyne demodulation scheme to get its in-phase and quadrature 
components J(t) and Q(t) respectively (for sake of simplicity details of the room 
temperature demodulation scheme, including several stages of amplification are 
not shown). Switching the signal ON and OFF with fms period allows to substract 
the noise background coming from the amplifier. 
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of this state. This continuous measurement does not induce spurious qubit excitation 
as long as the the intra-resonator photon number (n ~ 2.5 in our measurements) is 
much below the critical photon number n crit = {uo ge — u c ) 2 /Ag 2 ~ 30 above which the 
dispersive approximation fails. After reflection on the resonator, the signal is routed 
through a circulator to a cryogenic amplifier and is then measured by homodyne 
detection at room-temperature, yielding the two field quadratures J(t) and Q(t). 
The thermal fluctuations of the qubit state induce some phase noise on the reflected 
microwave signal, and thus some noise on each quadrature X(t) (X — I,Q). 

We start by computing the power spectrum of the qubit thermal fluctuations. 
Assuming that the bath consists of a bosonic Markovian bath at temperature T, as 
expected for the impedance of the electromagnetic environment, the qubit dynamics 
at thermal equilibrium is described by a simple rate equation [14] 

Pee = -Pgg = Pee + Fn th (l - 2p ee ) (1) 

where T is the qubit energy relaxation rate, and n t h{T) = (exp(fuv c /kT) — is the 
mean photon number at temperature T. This yields a steady-state population of the 
qubit excited state pf' = , "i' 1 , or z t h = — rnj: — after conversion into spin units 
z(t) = 2p ee (t) — 1. The corresponding noise power spectrum S z (ui) can be computed as 
the Fourier Transform of the two-time correlation function C z (r) = (z(r)z(O)} which 
is 

C z (t) = 4 exp (-r(l + 2n th )r) [l - p%] pf e . (2) 
After Fourier transform, we obtain 

Note that these expressions are only approximate because the transmon is not 
a genuine two-level system but an anharmonic resonator with an infinite number of 
excited states. The previous expressions are thus only valid in the limit where the 
population of these higher excited states is negligible, which in our case is true up to 
temperatures around 100 mK. 

We model the effect of the qubit state thermal fluctuations on the measuring 
signal quadratures X(t) by assuming that the field inside the resonator follows 
instantaneously the qubit state. Here this assumption is justified by the large 
bandwidth of the resonator, obtained by chosing a large resonator input capacitor. 
The quadratures are then simply expressed as X(t) = X + (AX/2)z(t) + where 
X is the average reflected signal for a qubit fully unpolarized, AX is the change in X 
when the qubit changes state, and £(t) is the total output noise of the amplifier. In 
the experiment we measure the sum of the noises on both quadratures 

Sv,on(u) = S/H + S q (uj) = St(u) + (AV/2) 2 S z {w), (4) 

where S%(w) is the output amplifier noise power spectrum and (AV/2) 2 = (AI/2) 2 + 
(AQ/2) 2 is the detector sensitivity. This quantity has the advantage of being 
insensitive to drifts of the phase between the local oscillator used in the demodulation 
and the measurement signal. It is worth noting that the mere presence of a continuous 
measurement of the qubit state has no effect on the dynamics of thermal fluctuations 
because this dynamics is fully incoherent and entirely goverened by Markovian rate 
equations [H| . The situation is very different when the qubit is continuously measured 
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Figure 2. (a) Noise spectra acquired for several temperatures T c (color solid 
lines) and lorentzian fits (dashed black lines), (b) Thermal population of the 
TLS as a function of temperature: comparison of the experimental data (black 
dots) with the theoretical prediction (red dashed curve), (c) Relaxation times as 
a function of temperature extracted from the widths of the loretzian spectra (blue 
dots) compared to the predictions of the model discussed in the text (red dashed 
curve) taking T _1 = T\ ^OmK = 226 ± 7ns, independently measured in a pulsed 
experiment at 20 mK. 

while being coherently driven, in which case the dynamics changes from diffusive Rabi 
oscillations to quantum jumps when the measurement strength is increased [16] . 

We measure the detector output noise spectrum Sy (w) for a series of temperatures 
T c . Each spectrum is measured after waiting 15 minutes thermalization time once the 
cryostat reaches T c . We also verify that the sample is well thermalized by acquiring 
two noise spectra for each T c , one upon warming up and the second upon cooling down, 
which are found to be nearly identical. Each spectrum is acquired by sampling I(t) 
and Q(t) with 100MHz sampling frequency. Each 1024-point set of the sampled signals 
is Fourier transformed and the amplitude of this transform is squared to obtain the 
noise spectra Si(uj) and Sq(lu), which are then corrected for the variations of setup 
gains in frequency and between the / and Q channels. The resulting spectra are 
summed to form Sv,on(u)- Each 2.5 ms the measurement microwave is turned OFF 
to measure the noise background of the amplifier Sv,off{w) = S^(u>) and subtract 
it from the signal. The resulting noise spectrum SV(w) = Sv,on{^) — Sv,off{u) is 
averaged typically 10 6 times. 

As shown in Fig 2a, the measured noise spectra SV(w) have a Lorentzian shape, 
with an amplitude rapidly increasing with temperature. The amplitude A and width 
T\ of each spectrum are fitted with a Lorentzian model ATi/(T\ +uj 2 ). According to 
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Eqs[3]andH the model predicts that i\ = T(l + 2n th ) and A = AV 2 [l - p*' e l ] pf e , 
The detector sensitivity AV 2 is experimentally calibrated in the following way : using 
exactly the same setup, we ensemble-average V 2 (t) = I 2 (t) + Q 2 (t) under saturation 
of the qubit g — > e transition with a second microwave source at frequency Lj ge /2iT. 
This yields AV/2 = 2.76 ± 0.14 mV. In this way we can directly extract from the fits 
the thermally excited state population p ee and the relaxation rates Ti as a function 
of the cryostat temperature T c . 

The fitted population (dots in Fig. [2b) agrees with the theoretical average 
population p*g (red dashed curve), calculated assuming two sources of radiation : 
the thermal field corresponding to the temperature of the cryostat coldest stage T c , 
with an average of nth{T c ) photons, and the thermal field radiated by the 30 dB 
attenuator thermalized at the still temperature Ts =600±100 mK, and attenuated 
(22 ± 0.5 dB) at 20 mK, contributing with n th {Ts) /10 2 2 photons. At the lowest T c , 
we find a thermally excited state population of 1 ± 0.5%, corresponding to an effective 
temperature of 55 mK. 

At T c = 20mK, the relaxation rate Tj" 1 deduced from the width of the Lorentzian 
noise spectrum (see Fig. Et) is found to be in excellent agreement with the qubit 
relaxation time Ti,20mif = 226 ± 7 ns, measured in a standard pulsed sequence. 
However, at higher T c , we observe that the fitted width decreases, which implies that 
the qubit energy relaxation time increases with temperature. This couterintuitive 
result disagrees with our model which predicts a relaxation rate T(n t h) = T(l + 2n t h) 
(see Eq. [J) increasing with temperature due to stimulated emission by the thermal 
field, yielding the red dashed curve in Fig. Gb (calculated with T — T^I^k)- 
This indicates that the qubit is not only coupled to its electromagnetic environment 
but also to another type of bath, causing some additional damping with a different 
temperature dependence. Additional support for this idea is that the measured 
relaxation time at 20 mK (226 ns) is significantly shorter than the expected damping 
time due to relaxation into the the external impedance at zero temperature (600 ns) , 
which indicates the existence of a non-radiative energy decay channel. We finally 
note that a similar increase of the relaxation time with temperature up to 150 mK was 
directly observed in a superconducting phase qubit, and attributed to non-equilibrium 
quasiparticles in the superconducting metal electrodes [17]; a similar scenario might 
explain our results. 

In conclusion, we have determined the thermal population of a superconducting 
qubit coupled to a resonator, even without single-shot detection capability, by studying 
the noise spectrum of its measuring signal. The population measured is in good 
agreement with the estimated thermal radiation reaching the sample. We observe 
however an increase in the relaxation time with temperature in contradiction with 
this model. This points to the existence of unknown non-radiative decay channels as 
observed in other qubit experiments [17] . 
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